Half-metallic semiDirac point generated by quantum confinement 

in Ti02/V02 nanostructures 



Victor Pardo^'^ and Warren E. Pickett^ 

^Department of Physics, University of California, Davis, CA 95616 
^ Departamento de Fisica Aplicada, Universidade de Santiago de Compostela, E-15782 Santiago de Compostela, Spain 

Multilayer V02/Ti02 nanostructures (d^ - interfaces with no polar discontinuity) are studied 
with first principles density functional methods including structural relaxation. Quantum confine- 
ment of the half metallic VO2 slab within insulating Ti02 produces an unexpected and unprecen- 
dented two-dimensional new state, with a semiDirac point Fermi surface: spinless charge carriers 
are effective-mass like along one principal axis, but are massless along the other. Effects of interface 
imperfection are addressed. 



VO2 is a magnetic oxide that undergoes a metal- 
insulator transition [l| upon lowering the temper- 
ature through 340 K, accompanied of a symmetry- 
breaking structural transition from the high temper- 
ature metallic rutile phase The insulating state 
takes place via a dimerization of the V-V chains Q . 
The origin of this metal- insulator transition is the fo- 
cus of much recent theoretical activity and remains 
uncertain. It could be due to the formation of a 
Peierls state 0, @| ) oi' it could be driven by correla- 
tions 0,13 J or more likely may have some mixed ori- 
gin Ti02 is isostructural (in one of its phases) 
and is a d*^ non-magnetic insulator that is very im- 
portant industrially and is well understood. 

The interface (IF) between a correlated insula- 
tor and a band insulator has been recognized as 
fertile ground for new behavior [lol 11|. and the 
LaTiOa/SrTiOa (LTO/STO) IF involving the Mott 
insulator LTO has attracted much of the theoretical 
study to date [13, [IE [I3| • For IFs between band in- 
sulators, LaAlOg/SrTiOajLAO/STO) has received 
a great deal of attention [11, [H, [13, 111, [11] . In both 
cases there is a polar discontinuity across the IF, 
and this aspect has been expected to dominate the 
resulting behavior, and lead to unexpected phenom- 
ena. 

The (001) V02/Ti02 IF has been studied by pho- 
toemission spectroscopy (PES) 20], which found the 
IF is insulating when the VO2 substrate is insulat- 
ing. PES also has uncovered [21] spectral weight 
transfer in V02/Ti02 thin films indicating strong 
correlation effects even for conducting VO2. Much 
of the focus on this nanostructurc has been on tun- 
ing the VO2 metal-insulator transition temperature 
[22| . because of its potential technological applica- 
tions. It has been found that a minimum thickness 
of 5 nm of VO2 is needed to sustain a metal- insulator 
transition, for thinner VO2 layers the transition no 
longer occurs 122] (the VO2 layer remains conduct- 
ing). The explanation is that the insulating state 



requires a collective structural dimerization along 
the rutile c-axis that is inhibited by confinement for 
thinner layers. Since the IF is not polar and the lat- 
tice mismatch is small (1 %), structural relaxation 
is not expected to be severe. Any unusual behavior 
of this multilayer (ML) will require different micro- 
scopic mechanisms than have been uncovered before. 

In this paper we present a theoretical study of 
the electronic behavior of the ML nanostructures 
(Ti02)„/(V02)m, denoted (n/m), looking in partic- 
ular at the evolution of the conduction and magnetic 
properties with VO2 layer thickness. Since VO2 is 
the component that is potentially conducting, we fo- 
cus on thin VO2 layers which will incorporate any 
consequences of quantum confinement. The prop- 
erties are found to depend strongly on layer thick- 
ness and the effects of quantum confinement at small 
thicknesses give rise to a new electronic state for cer- 
tain MLs. 
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FIG. 1: (Color online) Structure of the 5/3 Ti02/V02 
supercell corresponding to growth along the (001) axis, 
which is the metal chain direction of the rutile structure. 
VI, V2, V3 label the V ion sites beginning from the one 
nearest to Ti02. (Due to a symmetry with respect to 
the center of the VO2 slab, the V layers are V1-V2-V3- 
V3-V2-V1.) 



Starting from an average rutile structure of [n/ni) 
MLs, we performed volume and c/a optimization 
and an internal relaxation of the atomic positions of 
all the atoms. The main modification is lattice strain 
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FIG. 2: (Color online) Total density of states in the 5/3 
multilayer, showing the location of the V and Ti bands 
relative to those of O (which are unpolarized) . The V 
occupied majority spin bands (plotted upward) lie within 
a 2 eV gap in the minority spin. 



along the c axis, interpolating between the slightly 
different c lattice constants. We have studied sev- 
eral MLs, varying both Ti02 and VO2 layer thick- 
nesses, following identical procedures. Our elec- 
tronic structure calculations were performed within 
density functional theory 241 us ing the all-electron, 
full potential code wien2k |25j based on the aug- 
mented plane wave plus local orbital (APW-l-lo) ba- 
sis set [26| . The exchange-correlation potential uti- 
lized to deal with possible strong correlation effects 
was the LSDA-I-U scheme 271 128| including an on- 
site U and J (on-site Coulomb repulsion and ex- 
change strengths) for the Ti and V id states. The 
values U=3.4 eV, J= 0.7 eV have been used for both 
Ti and V to deal properly with correlations in this 
multilayered structure; these values are comparable 
(slightly smaller) than what have been used for bulk 
VO2 [1,0, [lOl. In the paper spin-orbit coupHng has 
been neglected. 

While we have studied a variety of n/m (001) 
MLs, we focus primarily on the (5/3) ML (1.5 nm 
Ti02, 0.9 nm of VO2). The structure, and the iden- 
tification of the three distinct V sites, is shown in 
Fig. [TJ In terms of distance from the Ti02 layer, 
the V ions are labeled VI, V2, VS. The tetragonal 
symmetry of the rutile structure has been retained 
in the x — y plane. 

V 3d bands (Fig. \^ dominate the spectrum close 
to the Fermi level i^p): and only three Ti02 cells 
are required to give negligible dispersion and thus 
confine the 3d states to a two dimensional (2D) sys- 
tem. FM alignment of the spins is preferred, and 



half metallicity results. Enlargement of the density 
of states (DOS) (not shown) reveals vanishing of the 
DOS precisely at Ei?, with no VI participation just 
below Ef. This curious vanishing of the DOS re- 
flects a zero gap semiconductor involving V2 and 
V3 ion states. 

The majority spin band structure along high 
symmetry directions shown in Fig. [3] clari- 
fies an unexpected and unprecendented electronic 
state. Two bands cross the Fermi level at a sin- 
gle point along the zone diagonal at the point 
fcs_D = (±0.37, ±0.37)7r/a (the precise position of kso 
along the (1,1) direction depends on the value of 
U). Inspection throughout the zone confirms that 
this Fermi surface crossing is a single point (rather, 
four symmetry related points), as is the Dirac point 
in graphene [3l| . This single point determines the 
Fermi energy, again as in graphene (3^ . The crossing 
of the bands precisely at E^? is therefore not acciden- 
tal, rather it is topologically determined: there are 
exactly six filled bands below this point, containing 
the majority spin electrons of each of the six V ions 
in the cell. 

These two bands crossing E^? involve separately 
V2 and V3 ion id states, as is illustrated with the 
color coded fatbands in Fig. [3] With no contribution 
from the IF ion VI, the dominance of the interior 
ions V2 and V3 identifies this as a quantum con- 
finement effect rather than an IF phenomenon. Ad- 
ditional VO2 layers, which relieve the confinement 
effects, add more bands and introduce a Fermi sur- 
face. In Fig. [3]we provide a surface plot of these two 
band energies in a small region in k-space centered 
on the band crossing point kgij. This state results 
only after the ion positions are relaxed, and arises 
due to band reordering at k—Q that occurs during 
the relaxation. The surface plot reveals yet another 
peculiarity: while the dispersion is linear along the 
(1,1) direction as is clear from the band plot, the 
dispersion is quadratic perpendicular to the diago- 
nal; the gap opens due to the loss of symmetry of 
the two bands off the diagonal k^ — ky, and does 
so quadratically. To differentiate this point from 
the graphene Dirac point, we refer to it as the (half 
metallic) semiDirac {sD) point. The corresponding 
mass tensor shows extreme anisotropy (zero to nor- 
mal values), as does the velocity (1.5 x 10^ cm/s 
to zero) [31|. This very strongly anisotropic disper- 
sion, between extremes (normal values, to zero) will 
give rise to peculiar transport and thermodynamic 
properties, which will be reported separately. 

The constant energy surfaces for both electrons 
and holes are plotted in the right panel of Fig. [3] 
for low energies in a small region around kso- The 
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FIG. 3: (Color online) Left panel: Blow-up of the band structure around the Fermi level showing in different colors 
the biggest character of each band. Notice that the two bands crossing at the Fermi level have character from the 
two most inner V atoms. Note the semiDirac point along the (110) direction where bands cross precisely at the 
Fermi level. Middle panel: Two different views of the same 'surface' plot of the two bands that cross the Fermi level, 
centered around the semiDirac point. The valence and conduction bands cross at a single point. The linear dispersion 
can be seen in the upper plot (upper left and lower right); the quadratic dispersion is clear in the lower panel, where 
the flatness of the conduction band is also clear. Right panel: contour plots at constant energy (in eV, relative to the 
Fermi level) of the valence band (below), and the flat conduction band (above) that leads to large M-centered Fermi 
surfaces for electron doping. 



conduction band has a flatness that opens a path 
for a Fermi surface to develop as a ring around the 
M=(7r,7r) point at very low electron doping. The 
valence band shows iso-energetic contours with an 
elliptical shape, with the longer axis perpendicular 
to the zone diagonal. 




FIG. 4: (Color online) Spin density plot, isosurface at 
0.15 e/A^. The particular orbital ordering in a spin- 
aligned configuration is shown. VI and V2 have one 
electron in a dy and V3 is in a dx orbital. 

In bulk VO2 (V^+i cations) the distortion from 
cubic symmetry of the VOe octahedron introduces 



a crystal field (actually, ligand field) that lifts the 
degeneracy of the t2g orbitals, splitting them into a 
d|| singlet and two d_L orbitals (using Goodenough's 
notation [4]). The orbital ordering that arises in this 
5 /3 ML is illustrated by the spin density isosurfaces 
in Fig. m The V ions (VI and V2) that terminate 
a V-V-V chain have an occupied dy orbital, whereas 
the chain-center V3 ion has an occupied orbital of dj^ 
of combined dxz and dyz character. This orbital or- 
dering is dependent on the magnetic ordering; when 
the spins along the c direction are antialigned (tit) 
the d|| orbital becomes occupied in all sites. 

Influence of VO2 slab thickness. A metal-insulator 
transition is observed [1^ for VO2 thicknesses above 
5 nm (approximately 15 layers), but experimental in- 
formation on crystalline samples with smaller thick- 
ness is sparse [33|. Our calculations show that the 
system has an insulating ground state for two layers 
of VO2, where spin antialignment is favored. How- 
ever, for a thickness of approximately one nm (3 
layers), the material is in the intermediate zero-gap 
semiDirac state described above, on the brink of 
metallicity. Thicker VO2 layers (four or more) be- 
come half metallic, a property that is much sought in 
oxide nanostructures because of its potentially enor- 
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mous technological applications in spintronic de- 
vices. 

Influence of magnetic alignment. We have stud- 
ied antialignment of the moments (ferrimagnctism) 
along the (001) V-V chains. Such AF coupling is 
energetically unfavorable in almost all cases. Inter- 
estingly, such antialignment changes the orbital or- 
dering: in the 5/3 multilayer (the semiDirac point 
system when FM), flipping the spin of the interme- 
diate V ion (V2) results in all V ions having an occu- 
pied d|| orbital, because the cr-bond along the z-axis 
between neighboring d|| orbitals favors AF coupling. 

Role of V-Ti exchange disorder. States that are 
very sensitive to disorder are less likely to have 
importance in applications, since thin film growth 
docs not result in perfectly ordered materials, so 
we have begun study of the effect of V/Ti ex- 
change near the IF. We find that the most un- 
expected feature, i.e. the development of a half 
metallic semiDirac point for three VO2 layers, is ro- 
bust with respect to two types of ion exchange that 
do not change the electron count, i.e. no doping. 
The first type was the interchange of VI with Ti 
across the IF, which is a typical defect in growth. 
If we label the Ti sites across the multilayer as 
Til /Ti2 /Ti3 /Ti4/Ti5 /Ti5 /Ti4 /TiS /Ti2 /Til , this 
first type of disorder corresponds to the interchange 
of VI and Til, corresponding to a non-abrupt IF. 
The second type is to interchange VI with Ti2, 
which are neighbors along the cation chain. In both 
cases a semiDirac point persists in spite of changes 
of the band structure, and confirms that it is the V2 
and V3 ions are produce the active bands. 

To address the robustness of this unusual property 
in the band structure for the 3-layer VO2 system, 
we have varied the thickness of the confining Ti02 
layer. Reduction of the Ti02 slab thickness to just 
three layers changes slightly the bands and thereby 
the position of the crossing point in the Brillouin 
zone, but still gives negligible dispersion along the 
z-axis, i.e. the behavior is still 2D. The semiDirac 
point is also robust with respect to the strength of 
correlation effects: the semiDirac point varies along 
the diagonal from (0.3,0.3) to (0.4,0.4) depending on 
both the choice of U on the V ions (for reasonable 
values, above 2 cV) and Ti02 thickness. 

The finding that quantum confinement, together 
with specific orbital occupation and perhaps impor- 
tant symmetries, in oxide multilayers can produce 
a semiDirac point at the crossover between insu- 
lating and conducting behavior introduces a novel 
feature in the physics of oxide hcterostructurcs: a 
polar discontinuity is not required to produce unex- 
pected and unprecedented electronic states in these 



systems. The transport behavior, and the changes 
with doping, for systems with a semiDirac point will 
be addressed in following papers, as will the compli- 
cating effects of spin-orbit coupling. We note that 
oxide nanostructures are mechanically more robust 
than graphene, which could make patterning of such 
multilayers possible. 
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